Analysis of genetic parameters and genetic trends for early growth and reproductive traits of Doyogena sheep managed under community-based breeding program

This study aims to estimate genetic parameters and genetic trends for early growth and reproductive traits of Doyogena sheep. Data used in the study were collected over 6 years (2013–2018). Studied traits were birth weight (BWT), weaning weight (WWT), 6-month weight (SMWT), average daily gains from birth to weaning (ADG0-3), average daily gains from weaning to 6-month age (ADG3-6), average daily gain from birth to 6-month age (ADG0-6), litter size (LS), lambing interval (LI), age at first lambing (AFL), and annual reproductive rate (ARR). (Co) variance components and genetic parameters were estimated using restricted maximum likelihood (REML). The analyses were carried out using WOMBAT program. Univariate analysis was applied to estimate genetic parameters. Six different animal models were fitted by including or excluding maternal effects. The direct heritability estimates for BWT, WWT, SMWT, ADG0-3, ADG3-6 and ADG0-6 were 0.33 ± 0.06, 0.31 ± 0.06, 0.14 ± 0.06, 0.13 ± 0.04, 0.11 ± 0.07, and 0.02 ± 0.05 respectively. Direct heritability for LS, LI, and AFL were 0.28 ± 0.12, 0.20 ± 0.5, and 0.001 ± 0.3, respectively. The maternal heritability estimates for BWT, WWT, and LS were 0.24 ± 0.12, 0.60 ± 0.07, and 0.24 ± 0.08, respectively. The genetic correlation between BWT with WWT and BWT with SMWT were 0.21 ± 0.07 and 0.21 ± 0.09, respectively. Genetic progress for most of the studied traits has shown promising improvements. Thus, continuation of selection, therefore, suggested for more improvements in the performance of Doyogena sheep. Direct heritability estimates decrease as lamb age increases and selection based on earlier body weight will be more efficient.


Introduction
Small ruminants, particularly native breed kinds play a significant role in the livelihoods of a considerable part of human population in the tropics from socio-economic aspects (Ahsani et al., 2010;Mohammadabadi 2021;Masoudzadeh et al., 2020a). Thus, combined trials with emphasis on administration and genetic progress to improve animal outputs are of decisive significance (Masoudzadeh et al., 2020b;Mohammadabadi et al., 2017). Economical and biological efficiency of sheep production enterprises generally improves by increasing productivity and reproductive performance of ewes (Mohammadabadi 2016;Amiri Roudbar et al., 2017, 2018Ghotbaldini et al., 2019;Mohammadabadi et al., 2021).
Community-based breeding programs (CBBPs) have recently attracted global interest as genetic improvement strategies in low input systems (Lamuno et al., 2018;Haile et al., 2019). CBBP has been designed to ensure the involvement of farmers (target groups) in all steps of the breeding program (Mueller et al., 2015). In Ethiopia, in 2012, the International Center for Agricultural Research in the Dry Areas (ICARDA), in partnership with the Southern Agricultural Research Institute (SARI), Areka Agricultural Research Centre (AARC) adopted CBBP in the Doyogena district to improve Doyogena sheep. Doyogena sheep are among the potential breeds of the country with better market preferences in the local market and Addis Ababa. The sheep have attractive morphological features with great potential for twining and fattening.
The potential genetic improvement of traits of interest was largely dependent on its heritability and genetic relationship among the traits of economic importance upon which selection may be applied. Information on heritability is essential for planning efficient breeding programs, and for prediction of response to selection (Falconer and Mackay, 1996). According to Berkana (2019), evaluation of any designed genetic improvement program is fundamental either to optimize the program if the designed improvement program is progressing towards the set goals or redesign other alternatives if it fails or deviates from the preset goals. Moreover, evaluation of genetic trends gives an indication of the genetic direction of the breed as well as the rate of genetic improvement from the time of application of the breeding program (Mallick et al., 2016). However, genetic studies of productive traits in sheep in low input systems are scarce due to lack of recorded data (Aguirre et al., 2017).
The present study has been planned to estimate the genetic parameters and trends for growth and reproductive traits and generate information for the optimization of the ongoing CBBP.

Description of the study area
The study area (Doyogena sheep CBBP) is found in Kembata Tembaro Zone of southern Ethiopia at a distance of 258 km to the Southwest of Addis Ababa (national capital). The district is located between 7 20 0 N latitude and 37 50' E longitude. Altitude ranges from 1900 to 2800 m. a.s.l. The average annual rainfall of the district is 1221 mm.

Breeding program description and animal management
Animals were identified by plastic ear tags. Enumerators were employed for routine animal identification, data recording, and followup. Enumerators use herd books for data recording. Selection of breeding rams takes place on a programmed date, twice a year. Researchers identify the candidate breeding rams by estimating best linear unbiased prediction (BLUP) breeding values for selection criteria traits using the performance and pedigree data recorded by the enumerators. Lambs wean at the age of 90 days (three-months age). Then the candidate rams pre-selection and ranking take place based on weaning weight. In the second stage-breeding ram ranking was carried out based on SMWT estimated breeding values (EBVs). Top (10% of the candidates) breeding rams were retained for breeding to be used in the community flock while the next best (positive EBVs) were sold for breeding purposes to other communities.
Culled males (negative EBVs) were either castrated or marketed to prevent unwanted mating (Haile et al., 2019). Selected best breeding rams usually serve not more than one year in the community flock. After one year of service, the breeding rams were sold to another area of the region.
The main feed sources for animals included Enset or false banana (Ensete ventricosum (Welw. Chessman) products of Amicho (fresh parts of false banana could be cooked like potatoes), corm (swollen underground parts of false banana), crop residue, improved forage/grass, kitchen leftover, and purchased concentrates. Flocks graze with tethering on the small private land. Free veterinary service was provided for CBBP participant farmers by ICARDA and SARI.

Data sets
The empirical data for the study were obtained from five ongoing breeder cooperatives. The performance data along with pedigree information were maintained in the data-recording book of individual breeder cooperatives. The data routinely collected by the enumerators were recorded at the time of the event. BWT was recorded within 24 h of lambing; WWT was taken at 90 days of age, and 6-month weight was taken at the age of 180 days. The average daily body weight gains from birth to weaning age; weaning to six months of age, and from birth to six months of age have been estimated as under: Average daily body weight ðBWÞgain up to weaning age Average daily BW gain from weaning up to 6 months age Average daily BW from birth to 6 month age Where: BWT ¼ Birth weight, WWT ¼ Three month weight/weaning weight at 90 days, SMWT ¼ 6-month weight at 180 days The annual reproductive rate, which is the number of lambs born per breeding ewes per year has been estimated according to Shigdaf (2013) as under: Where: ALS ¼ Average litter size; LI ¼ Lambing interval; 365 ¼ Days of one year;

Statistical analysis
The variance components and resulting genetic parameters were estimated on a model-fitting effect of parity, year of birth/year of lambing, season of birth/season of lambing, type of birth, sex, and sites (cooperatives) as fixed factors. Based on the first performed analysis of variance using SAS (SAS, 2009), significant fixed effects were identified to be included in the models. Then, significant fixed effects were fitted in the subsequent models for estimating genetic parameters. Genetic parameters for reproductive traits were estimated by WOMBAT software (Meyer, 2012). Univariate animal models were fitted to estimate genetic parameters. Similarly, bivariate, trivariate, and multivariate analyses were applied for estimation of correlations. Direct additive and maternal additive genetic effects with or without a covariance between them, and maternal permanent environmental effects were tested for all traits in different combinations to yield six models. The six models were as follows: (Co) variance components models y ¼ xb þ z1a þ z2m þ e; with covða; mÞ ¼ Aσam (model 4) y ¼ xb þ z1a þ z2m þ z3c þ e; with covða; mÞ ¼ 0 (model 5) y ¼ xb þ z1a þ z2m þ z3c þ e; with covða; mÞ ¼ Aσam (model 6) Where: y ¼ vector of observed traits; b, a, m, c ¼ Vectors of fixed effects, direct additive genetic effects, Maternal additive genetic effects, maternal permanent environmental effects and Vector respectively. e ¼ residual effects X, Z, 1 , Z 2, and Z 3 ¼ Incidence matrices, respectively relating fixed effects, direct additive genetic effects, maternal additive genetic effects, and maternal permanent effects to y; All components with the phenotypic variance (σ2p) being the sum of σ2a, σ2m, σam, σ2c, and σ2e, were derived at convergence. Depending on the model, we computed Direct heritability as; h 2 ¼ σ 2 a σ 2 p (11) Maternal heritability; m 2 ¼ σ 2 m σ 2 p and; (12) the direct-maternal covariance as proportion of phenotypic variance; (c am ¼ σ am /σ 2 p), with a corresponding estimate of the direct-maternal correlation [r am ¼ c am /(σ 2 a *σ 2 m )]. Similarly maternal environmental variance ratio was estimated by the maternal permanent environmental variance as a proportion of σ 2 p (p 2 ¼ σ 2 c/σ 2 p). The genetic correlation between direct and maternal genetic effects (r am ) is estimated as the ratio of the estimates of the σ am to the product of the square roots of the estimates of σ 2 a and σ 2 m (Meyer (2007).
Total heritability (h 2 t) was calculated according to the following equation (Willham, 1972).
To determine the most appropriate model, likelihood ratio tests (logL) were used for each trait. An effect was considered to have significant influence, when its inclusion caused significant increase in loglikelihood, compared to the model in which it was ignored. When loglikelihoods did not differ significantly (P > 0.05), the model that has fewer parameters was selected as the most appropriate model. All models included direct additive genetic effect and this was the only random factor in Model 1. Model 2 included the maternal permanent environmental effect, fitted as an additional random effect. Model 3 included an additive maternal effect fitted as second random effect. Model 4 was the same as Model 3 but allowed for a direct maternal covariance Cov (a,m) . Model 5 and Model 6 included additive maternal and maternal permanent environmental effects, ignoring and fitting, respectively, directmaternal covariance.
The genetic trends were estimated by the weighted regression of the average breeding value of the animals on the year of birth or year of lambing. These procedures were carried out with statistical program R core software (R, 2018) and MS-excel. Genetic change for the traits over the selection period was calculated by subtracting the mean of the estimated breeding values at the beginning of the CBBP from the mean of the EBV at the time of this study's 2018-selection year.

Results and discussion
3.1. Genetic parameter estimation 3.1.1. Data and pedigree structure In Table 1 the structure of pedigree has been presented.

Inbreeding
The coefficient of inbreeding (F) showed an increasing trend within the 6-year selection period ( Figure 1). The coefficient of inbreeding was assumed to be zero until the year 2014; afterward, it increased with the selection years. At the time of this study (2018 selection year) coefficient of inbreeding was 0.30% with an average annual inbreeding trend of 0.08%. The proportions of inbred animals from 4497-studied animals were 37 animals. The most likely reason for this inbreeding increment could be the selection of superior breeding rams without seeing their detailed pedigree. FAO (2010) suggested that the inbreeding rate must be maintained lower the range of 0.5-1% per year to avoid risks of genetic disorders and inbreeding depression. The inbreeding coefficient obtained for Doyogena sheep was considered as an acceptable percentage, however, F is in increasing trend, and thus consideration should be given during allocation of breeding rams. Negussie et al. (2002) reported 19 years (1978)(1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997) inbreeding coefficient of 0.78% with an annual trend of 0.07% for Horro sheep on the station management system. Gizaw et al. (2013) reported F for Menz sheep after 10-years selection was 1.7% with 0.17% increment per generation, which is higher than the presently estimated inbreeding coefficient. Generally, inbreeding leads to a reduction in additive genetic variance and heritability (Falconer and Mackay, 1996).

Genetic parameter estimates for growth traits
Based on the appropriate models the estimates of direct heritability (h 2 a) for BWT, WWT, and SMWT were 0.33 AE 0.06, 0.31 AE 0.06, and 0.14 AE 0.06, respectively (Table 2). Except for the moderate heritability estimate for SMWT, which reflects less genetic variation among lambs, the estimates of h 2 a for BWT and WWT fall within the range of values reported in the high heritability value. The present estimate of direct heritability for BWT (0.33 AE 0.06) was found in the range reported by Yacob (2008)    The estimate of h 2 a for WWT (0.31 AE 0.06) was also found in the estimated range reported by Yacob (2008) for Afar sheep (0.11-0.37) and BHS sheep (0.00-0.29) but lower than the estimate for Menz sheep (0.46) by Gizaw et al. (2007). The estimate h 2 a for WWT by Abegaz et al. (2002) for Horro sheep (0.16 AE 0.05) and Gizaw et al. (2014b) for Menz sheep (0.19) were lower than the present estimate.
The estimate of h 2 a for SMWT 0.14 AE 0.06 was found in the range of h 2 a estimated for Afar sheep (0.11-0.37) and BHS (0-0.29), while, the report of Gizaw et al. (2014b) for Menz sheep (0.46) were much higher than the present estimate. Abegaz et al. (2002) estimated 0.18 AE 0.05 of h 2 a for Horro sheep, which is lower than the current estimate. From the genetic point of view, h 2 a estimated for BWT and WWT indicated that high variation within the breed would be a greater opportunity for selection response during genetic improvement through selection for these traits. Moreover, WWT will be the best criterion for selection to increase the pre-weaning growth rate because, selection based on BWT, which has the highest heritability could cause dystocia. However, the confounding effect of direct genetic and maternal genetic effects needs care.
The permanent maternal environmental effect (c 2 ) for BWT was moderate in this study (0.20 AE 0.09). This indicates the importance of maternal environment and care at the birth of lambs. The current estimate is similar to the findings of Gowane et al. (2010) in Bharat Merino sheep (0.19 AE 0.02). For SMWT trait, the maternal environmental effect is more important than maternal genetic effect. Venkataramanan (2013) reported a similar finding for Nilagiri and Sandyno Indian sheep breeds. The result suggested that, even if maternal effects tend to diminish with age, some adult traits would nevertheless contain this source of variation. The current finding indicated that maternal heritability (h 2 m) is an important ratio for BWT, WWT, and the estimates were 0.24 AE 0.12 and 0.60 AE 0.07 respectively.
In the current study, for BWT, h 2 m estimate was in the range estimated (0.06-0.46) for blackhead somilan sheep (BHS) by Yacob (2008). Abegaz (2002); Assan et al. (2002) and Yacob (2008) estimated h 2 m of 0.12 AE 0.2, 0.24, and 0.02-0.21 for Horro sheep, Sabi sheep, and Afar sheep respectively and these all value are lower than the present estimate. However, Mousa et al. (2013) reported higher estimated values of h 2 m for Moroccoian Timahdit sheep (0.59). The present estimated  maternal heritability of 0.6 AE 0.07 for WWT was higher than the above-mentioned maternal heritability estimates. High and negative additive maternal genetic correlation estimates were observed (Table 3) for BWT (À0.61 AE 0.15) and WWT (À0.81 AE 0.11) traits. Similar results were summarized by Safari and Fogarty (2003) for a wide range of sheep breeds. The correlation estimates between direct additive and maternal genetic effect (r am ) for both the traits become negative means improvement in one will result in a reduction of another. The result might be due to the structure of the data set used in the analysis i.e. the number of generations the animals were measured both directly and as dams were limited causing lack of large pedigree. The estimates of total heritability (h 2 t) for BWT, WWT, and SMWT were 0.21, 0.12, and 0.14, respectively. The h 2 t estimated by Abegaz (2002) for Horro sheep for BWT, WWT, and SMWT were 0.14, 0.12, and 0.21, respectively showing little increment across lamb age, which is slightly in contrast with the present results. The result indicated that maternal effects were important for weights until about 6 months of age. Based on the best-fitted models, the estimated h 2 a for ADG0-3, ADG3-6 and ADG0-6 were 0.12 AE 0.04, 0.11 AE 0.07, and 0.02 AE 0.05 respectively ( Table 3). The estimate indicated that the inclusion of maternal permanent environmental effects in the analyses could improve the models for daily weight gain traits. The fractions of maternal permanent environmental variance are highly reflected for all considered average daily weight gain traits.
The estimate indicates variance due to permanent maternal environmental effects (c 2 ) for ADG0-3 (0.21 AE 0.03) and ADG0-6 (0.26 AE 0.04) have been found significantly higher than later age daily weight gain traits of ADG3-6 (0.09 AE 0.04). It decreases with increasing lamb age. This could be due to the influences of feeding levels at the later age of the lambs and the maternal behavior of the dam, especially for preweaning growth traits in the lambs. The value of maternal permanent environmental variance in model (2) for this trait is not significantly different from other models' values. The result is also found in the range reported by Yacob (2008) for Afar and BHS sheep and lower than the report of Radwan and Shalaby (2017) and Matika (2001) for Rahmani and Sabi sheep, respectively. The estimate of h 2 a for ADG3-6 was comparable with the report of Yacob (2008) which is 0.00 and 0.09 for Afar and BHS sheep under station management conditions. σ 2 a ¼ direct additive genetic variance; σ 2 c ¼ maternal permanent environmental variance;σ 2 m ¼ maternal additive genetic variance; σ am ¼ additive and maternal additive genetic covariance, σ 2 e ¼ residual variance, σ 2 p ¼ phenotypic variance, h 2 a ¼ direct heritability c 2 a ¼ ratio maternal permanent environmental variance to phenotypic variance, h 2 m ¼ maternal heritability; r am ¼ correlation between direct maternal additive genetic effects, h 2 t ¼ total heritability and SE ¼ standard error.

K. Habtegiorgis et al. Heliyon 8 (2022) e09749
The estimates of h 2 t values for ADG0-3, ADG3-6, and ADG0-6 were 0.12, 0.11, and 0.023 respectively, which is in a similar range to the h 2 a estimates. The estimates were in the moderate range except for ADG0-6. Total heritability estimates for ADG0-3 and ADG0-6 were comparable with the finding of Abegaz (2002) for Horro sheep, which were 0.13 AE 0.04 and 0.04 AE 0.03, respectively.

Genetic parameter estimates for reproductive traits
The estimate of h 2 a and h 2 t for litter size were 0.28 AE 0.12 and 0.29, respectively (Table 4). The h 2 a estimate indicated that genetic improvement through direct selection for this trait would be high for Doyogena sheep. Compared with another study the current heritability estimate for litter size was higher. Abegaz (2002) estimated for litter size using direct additive and repeatability models were 0.15 and 0.07 respectively for Horro sheep. Matika (2003) and Mohammadi et al. (2012) reported estimated heritability of 0.26 for Sabi and 0.14 for Zandi sheep by fitting the threshold model. Haile et al. (2020) reported h 2 a of 0.08 AE 0.041 and 0.08 AE 0.04 for Bonga and Horro sheep breeds respectively under CBBP and that is lower than the current estimates. Likewise, Khan et al. (2017), and Mohammadabadi and Sattayimokhtari (2013) estimated litter size heritability of 0.25 and 0.06 AE 0.02 for Rambouillet and Kermani sheep breeds in India and Iran respectively and the h 2 values reported were lower than the current finding.
The h 2 t of 0.29 means that 29% of the variability in this trait between the ewes is due to genetic differences among ewes. Therefore, the current estimate of heritability shown for litter size has the highest heritability σ 2 a ¼ direct additive genetic variance; σ 2 c ¼ maternal permanent environmental variance; σ 2 m ¼ maternal additive genetic variance; σ am ¼ additive and maternal additive genetic covariance, σ2e ¼ residual variance, σ 2 p ¼ phenotypic variance, h 2 a ¼ direct heritability c 2 a ¼ ratio of maternal permanent environmental variance to phenotypic variance, h 2 m ¼ maternal heritability; r am ¼ correlation between direct maternal additive genetic effects, h 2 t ¼ total heritability and SE ¼ standard error. Under the best-fitted model, the estimated h 2 a and h 2 t for LI were 0.20 AE 0.5 and 0.20 respectively. Lobo et al. (2009) and Abdoli et al. (2019) estimated h 2 a of 0.06 and 0.02 for LI for Brazilian multi-breed sheep and Iranian Lori-Bakhtiari sheep and the values reported were lower than the current estimate. The h 2 a and h 2 t estimate for AFL were 0.001 AE 0.32 and 0.007, respectively. Since AFL is strongly influenced by environmental effects, low heritability estimate was obtained. Selection based on AFL performance may result in slow genetic improvement. Therefore, selection for ewes for AFL trait of ewes should be based on female relatives of ewes or on correlated traits, which have high and positive genetic correlation with ewes AFL. Lobo et al. (2009) and Abdoli et al. (2019) estimated 0.04 and 0.07 of AFL h 2 a for Brazilian multibreed meat sheep and Lori-Bakhtiari sheep breeds respectively. The h 2 values reported were higher than the current finding.

Correlation estimates between reproductive traits
The bivariate analysis of genetic correlation between reproductive traits is given in Table 5. Negative estimate of genetic correlation -0.44 AE 0.9 was obtained between LI and LS. However, the genetic correlation between AFL and LS was strong and negative (À0.98 AE 0.32). The present estimate of genetic correlations was similar to the other estimates reviewed by Safari et al. (2005). Khan (2017) also reported comparable genetic correlation of -0.00 AE 0.02 and -0.006 AE 0.02 for LS with AFL and LS with LI in Rambouillet Sheep respectively.

Genetic correlations between growth traits
The present study indicated that BWT had weak genetic correlation with the studied body weights and daily weight gain traits ( Table 6). The negative correlation (À0.35 AE 0.14) between AD'G0-3 and ADG3-6 indicated that lambs that grew faster in the pre-weaning period, grew more slowly during post-weaning period and vice versa. Mohammadi et al. (2015) reported similar findings to the present study in the Lori sheep breed. A moderate positive genetic correlation was observed between WWT and SMWT (0.52 AE 0.09). The positive genetic correlations between the two traits indicated that the genes that are responsible for increasing WWT result in increasing in SMWT trait. It could be used as selection criteria for improvement in body weight traits. The positive and moderate genetic correlation that post-weaning body weights and body weight gains may be under the influence of the same set of genes (Pleiotropy). In the ongoing CBBP both WWT and SMWT traits were considered the most appropriate selection criteria (Jembere et al., 2016). The current genetic correlations were similar to the report of Abegaz (2002) estimated for Horro sheep.

Birth weight (BWT)
Both direct genetic and maternal genetic trends had fluctuating trends ( Figure 2). The annual estimated direct genetic trend was negative (À0.0026 kg/year) and insignificant (p > 0.05), however, the maternal genetic trends showed an increasing trend (0.0023 kg/year. The estimated direct and maternal genetic changes were 0.00085 kg and -0.004 kg respectively. When compared to other studies, Gizaw et al. (2014a) reported higher and positive genetic change in BWT for Menz sheep. The author reported a genetic gain of 0.005 kg in the 4th generation of selection program. Gholizadeh et al. (2015) reported an annual genetic gain of 0.00 kg/year for Baluchi sheep, which is almost similar to the current finding. Since direct genetic gain for the traits showed slightly negative trend, demonstrates that these traits should not be taken into consideration in the selection process by breeder cooperatives. Figure 3 shows the value of direct genetic and maternal genetic trends over the selection period. The magnitude of direct genetic trends Table 6. Estimates of genetic (below diagonal) and phenotypic (above diagonal) correlations between growth traits in multi-trait analysis.  estimated illustrates that there has been significant (p < 0.05) and positive genetic improvement in WWT with 0.3 kg in a period of 6-years selection (0.048 kg/year). From 2015 to 2017, the direct genetic trend was in decreasing trend, after which there was a slightly increased in values. The negative trend obtained could be due to the reason that less rainfall caused lack of forage and then decline in the growth performance of the flock within the selection years. The maternal genetic trend had a decreasing trend by -0.051 kg/year and -0.23 kg per 6-years selection. The direct genetic trend of 0.048 kg/ year was higher than the study of Mostafa et al. (2011) for the Arman sheep breed (0.007 kg/year) and lower than the report of Mokhtari (2010) for kermani sheep (0.125 kg/year). The direct genetic change (0.3 kg) was positive and significant for WWT. Gizaw et al.(2014a) reported 0.45 kg of genetic change in the 4 th generation of Menz sheep which is higher than the current finding.

Weaning weight (WWT)
In Figure 4 and Table, 7 annual genetic trend and genetic change for SMWT has been shown. The existing method of breeding ram selection was based on SMWT. The estimated annual direct genetic trend (0.036 kg/year) was positive and highly significant (p < 0.01). The fit of the regression shows 73.4% coefficient of determination with the regressed value. Direct genetic change was 0.15 kg. The estimate of direct genetic

Breeding value (kg)
Birth year WWT-maternal genetic trend    change (0.151 kg) for SMWT provides a good picture of the selection program concerning SMWT and therefore, continuation of selection based on SMWT trait is suggested to become more successful. The present estimate of direct genetic trend was concurrent with the study of Mokhtari and Rashidi (2010) for Kermani sheep (Mohammadi et al., 2011) and for Zandi sheep was 0.021 kg/year. Shaat (2004) and Arora et al. (2010) reported higher estimates in Rahmani sheep (0.135 kg/year) and for Malpura sheep (0.061 kg/year respectively. The present estimate is lower than the report of Gizaw et al. (2014a) in Menz sheep; they report 1.3 kg genetic progress at the 4 th generation.

Reproductive traits genetic trends
Direct genetic trend for litter size was fluctuations ( Figure 5) and nonsignificant (p > 0.05). This fluctuated and slow genetic trend, for litter size, could be because of lack of data quality. Genetic progress for LS showed improvement across selection years and could be taken into consideration in the process of selection.
LI shows ( Figure 6) a decreasing trend (À0.0007 days/year). However, AFL trait showed an increasing trend (0.0174 days/year). Annual genetic changes in each of the reproductive traits considered in the present study are presented in Table 8.

Conclusion
The objective of this study was to evaluate community-based breeding program breed improvement strategy to further improve genetic gains in Doyogena sheep. Promising results of selection were observed from the ongoing Doyogena sheep CBBP. The different estimates of heritability obtained from the different models suggest that model choice is an important aspect of obtaining reliable parameter estimates to be used in genetic parameter estimation. The moderate to high estimated heritability for the growth traits suggested the scope for further improvement of these traits. Coefficient of inbreeding showed an increasing trend across selection years. There was a negative genetic trend observed for BWT trait. However, WWT and SMWT traits were genetically improved across the years. The estimate of genetic gain for SMWT trait was the greatest among the body weight traits. The result suggested that selection for litter size trait would be more efficient, due to its higher genetic variance; however, it may hurt growth performance.

Author contribution statement
Kebede Habtegiorgis: Conceived and designed the experiments; Analyzed and interpreted the data; Wrote the paper. Aynalem Haile; Tesfaye Getachew: Performed the experiments; Contributed reagents, materials, analysis tools or data; Wrote the paper.
Manzoor Ahmed Kirmani: Performed the experiments; Contributed reagents, materials, analysis tools or data.

Funding statement
This work was supported by the International Centre for Agricultural Research in the Dry Areas (ICARDA) and Southern Agricultural Research Institute (SARI).

Data availability statement
Data will be made available on request.

Declaration of interest's statement
The authors declare no conflict of interest.

Additional information
No additional information is available for this paper. national Center for Agricultural Research in the Dry Areas (ICARDA) for allowing data and providing technical support for the researcher. We are also grateful to all persons who were involved in this activity for their contribution in different aspects.